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Most of the reflectometry methods which are used for determining the phase of complex reflection coefficient 
such as reference method and Variation of Surroundings medium are based on solving the Schrodinger equa- 
tion using a discontinuous and step-like scattering optical potential. However, during the deposition process 
for making a real sample the two adjacent layers are mixed together and the interface would not be discon- 
tinuous and sharp. The smearing of adjacent layers at the interface (smoothness of interface), would affect the 
reflectivity, phase of reflection coefficient and reconstruction of the scattering length density (SLD) of the sam- 
ple. In this paper, we have investigated the stability of reference method in the presence of smooth interfaces. 
The smoothness of interfaces is considered by using a continuous function scattering potential. We have also 
proposed a method to achieve the most reliable output result while retrieving the SLD of the sample. 
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1. Introduction 

In the past decades, neutron reflectometry has been developed as an atomic-scale probe with appli- 
cations to the study of surface structure ranging from liquid surfaces to the solid thin films. The type and 
thickness of the nanostructure materials can be determined from measuring the number of neutrons, 
reflected elastically and specularly from the unknown sample HjiSJ. 

Neutron reflectometry problems are generally divided into two groups; "Direct problems" and "In- 
verse problems". The goal of the direct problems is to solve the Schrodinger equation for a distinct sample 
so that the neutron wave function is determined; On the other hand, the application of inverse problems 
is to extract the information of the interacting potential by using the complex reflection coefficient 1 1 1. 

Measuring the intensity of reflected neutrons in terms of the perpendicular component of neutron 
wave vector to the surface, q, provides us with useful information about the scattering length density of 
the sample along its depth. As with any scattering technique in which only the intensities are measured, 
the loss of the phase of reflection would result in unresolvable ambiguities on retrieving the SLD from 
measurement. Without the knowledge of the phase of reflection, more than one SLD can be found for the 
same reflectivity data. By knowing the phase of the reflection, a unique result would be obtained for the 
depth profile of the unknown sample |4 5 1. 

In the recent years, several methods have been worked out for determining the complex reflection co- 
efficient such as: Dwell time method, variation of surroundings medium and the reference layer method 
which seems to be the most practical |1 1. 

All of these methods are worked out for an ideal sample in which the interface of two adjacent layers 
is discontinuous and sharp. In this case, the reflection coefficient is determined by the solution of one 
dimensional Schrodinger equation for a step-like optical potential. However, as we know from a real 

* E-mail: s.jahromi@dena.kntu.ac.ir 



© S.S.Jahromi, S.F. Masoudi, 2012 13604-1 



S.S. Jahromi, S.F. Masoudi 




(a) (b) 



Figure 1. Smearing of the SLD of two adjacent layers at the interface. The SLD of layer have been mixed 
together during the deposition process, (a) Smearing of layers in atomic scale, (b) Continuous variation 
of SLD at the interface in one dimension. 



sample, there is some smearing at the boundaries. This generally happens during the deposition of a 
top layer which is miscible with the bottom material (figure [ij. This process is strongly temperature 
dependent. In this case, the interface is defined as a thin layer across which, the SLD of the sample varies 
smoothly around the mean position of the interface | 6 |. 

In most of the simulation methods, the reflectivity measurement and phase determination process 
is performed for an ideal sample. Considering a continuous and smooth varying optical potential at the 
interface, it would affect the retrieved reflection coefficient. In the present work, we have studied the ef- 
fects of potential smoothness on retrieving the depth profile of the sample and the phase of the reflection 
by considering a continuous potential at boundaries, using reference layers method. As retrieving the 
SLD from the real and imaginary part of the reflection coefficient is very sensitive, we are going to inves- 
tigate the way one can find the characteristics of an unknown sample in the presence of the smoothness 
effects. 



2. Theory and method 

Interaction of neutron with nuclei is demonstrated by the neutron optical potential viz), which is a 
function of the scattering length density of the sample, v{z) = 2nh^p{z)/4n where p{z) is the scattering 
length density (SLD) profile as a function of the coordinate z; normal to the sample surface EE). As 
the thickness of the sample is of the order of nanometer, multiple scattering is neglected in neutron 
specular reflectometry and the neutron specular reflection is accurately described by a one-dimensional 
Schrodinger wave equation 

{dl +[q^- Anp {z)\}^{q,z) = 0, (1) 

where y/iq,z) is the neutron wave function and qo is the z-component of the incident neutron wave 
vector in vacuum. 

Specular reflection is rigorously invertible, however, if both the modulus and the phase of the com- 
plex reflection coefficient, riqo), were known, a unique result for scattering length density, piz) would 
be determined |5 |. 

An exact relation for the reflection coefficient riqo) can be derived from the transfer matrix which 
is a 2 X 2 matrix with the elements A{qo),B{qo),C{qo), and D{qo) which are determined by solving the 
equation for a film with SLD fSllT). 

Correspondingly, for a sample with arbitrary surroundings, the solution of equation jlj in terms of 
the transfer matrix elements is represented by: 

where tiq) is the transmission coefficient and / and h are the refractive index of the fronting and backing 

medium having SLD values of pf and p^ respectively. For non-vacuum surroundings, n=\l- ^n^j 

where n stands for / and h. Correspondingly, for vacuum fronting or backing niq) = f = h = l 
The transfer matrix has a unit determinant and is unimodular; 

Aiqo)Diqo)-Biqo)Ciqo) = l. (3) 
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As the other property of the transfer matrix, reversing the potential, piz) — ^ p{L- z), would cause the 
interchange of the diagonal elements A and D without having any effect on the off-diagonal elements, B 
and C 111. 

By solving the equation for the reflection coefficient, r{q), we have: 

_ [fh^B^ + fD^) - [h^A^ + C2) - 2i [fh^AB + fCD) 
[f^h^B^ + /2D2) + [h^A^ + C2) + 2fh * 

The reflectivity Riq) = \ riq)\^ can be represented in terms of the transfer matrix elements by introducing 
a new quantity I.{q): 

Uq) = \^^fh = [fh^B^ + f^^) + (^^^^ + (5) 

where R is the amplitude of the complex reflection coefficient and the superscript fh, denotes a sam- 
ple with fronting and backing medium with refractive index / and h, respectively. By introducing the 
following new quantities 

af^ = [r'hA^^f-'h-'C^]^ 
l3f^ = [fhB^ + fh-^D^], 

rf^ = hAB-^h~^CD. (6) 
Equations l|4| and l[5j, are represented as 

rJ^ = — ^— , (7) 

Ziq) = l3f^ + af^. (8) 

These two equations are the two key relations in reference method. Equation ^ shows that below the 
critical wave number qc = 47TMax[gf,ghY^^, the method is not applicable since the information of riq) 
is not available. However, the missing data below the qc could be accurately retrieved by extrapola- 
tion CIItI. In section 4 we will introduce an algebraic method for determining the missing value below 
the critical q. 

2.1. Reference layer method 

Here we explain how one can find the complex reflection coefficient of any unknown non-absorptive 
layer by using some reference layers. Suppose a known layer has been mediated between the fronting 
medium and the unknown layer. In this case, the transfer matrix of the whole sample (containing the 
known and unknown layers) can be expressed by the multiplication of the transfer matrix of each one as 
follows: 

(A B\_(a b\(w X 
[c D)-[c d)[y 

where (ay... yd) refer to the matrix elements of the unknown part and (w;, . . . , z) present the matrix ele- 
ments of the known part. Consequently we have 

Z = [h^a^ + c^] [w^ + fx^] -h [h^b^ + d^] [y^ + fz^] + 2 [h^ah + cd) [wy + fxz] . (10) 

Alternatively, the above formula can be denoted in compact notation in terms of three unknown real- 
valued parameters as 

Z(^) = 2i±|^ = h^^a'^'^ + fai^^u + 2//zf^^rS'' > dD 

where the subscripts K and U refer to the known and unknown parts of the sample and the tilde rep- 
resents the parameters of the sample with mirror reversed potential {A ^ — > D). ff and hh denote the 
same fronting and backing media having the SLD values of pf and p^, respectively. 



(9) 
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In order to determine the reflection coeflicient, riq), we consider three measurements of reflectiv- 
ity for the three known reference layers. Using equation j6) and (TT) for three different reference lay- 
ers (Ki,K2,K3), three equations will be obtained. By solving these equations in matrix algebra as shown 



in equation |12|, the parameters (a^ , ,7^ ) are determined which represent the r"'"'iq) for the un- 
known film in contact with a uniform identical surroundings in front and back with p^- 



M: 



f2^ff 
J "K3 



Ki 



2fhf/V' 



(12) 



For q> qc this gives us the reflection coefficient 

fihh _ ^hh _ n\n.rhh 



(13) 



which is equivalent to the reflectivity of a "free fllm" deflned by PequivC-z) = p^iz) - p^-By inverting 
and shifting by a known constant value, the SLD of the free fllm can be retrieved 1 1 8 |. 



hh 



2.2. Determining the phase below the critical wave number 

There are several ways for retrieving the missing data below the critical wave number such as ex- 
trapolation or using polarized neutron and a magnetic substrate. Here we introduce an algebraic method 
to retrieve the data below qc. 

The transfer matrix for a layer having constant SLD in terms of wave number is represented by: 



cos nqoL j^sinnqoL 
-nsinnqoL cos nq^L j' 



(14) 



where n = yjl-^ngl q^, L is the thickness and g is the SLD of the sample. 

Introducing several new parameters, peq = Q-Qh^ ^eq = \J ^^Qh ^ricl ^eq = yjl-4ngeq/ ql^^, for 



T we have 



(15) 



cos rieq qeqL sin rieq qeqL 

[-hrieq sin Ueq qeqL COS Ueq qeqL 

By inserting the elements of transfer matrix in equation ( [13) and using equation equation l[7|, we 



have: 



(16) 



Equation 1I61 shows that knowing r^^{g) for qo ^ y/4ngh, the reflection coefficient of a free sample 

having SLD value of ^eq = Q-Qh^ r^^io - Qh)^ known for specific wave number qj = \J q^- ^^gh for 
qj^O. 



2.3. Smooth variation of interface 



As we mentioned in section [2TT] there is some smearing at boundaries in real samples and the SLD 
varies smoothly and continuously from one layer to another. The smooth variation of SLD at boundaries 
can be expressed by Error function as follows |9 1: 



g{z) = gi{z) + 



g2iz)-gi{z) 



1-herf 



z-A\ 
V2al 



(17) 



where a is the smoothness factor. In other word, a is the thickness of the smoothly varying area across 
the interface. A is also the mean position of the interface or the turning point of the Error function. Since 
the reference layer method is based on discontinuity of the interfacial potential, it is obvious that the 
effect of smoothness at boundaries would affect the reconstruction procedure. In the next section we 
have investigated these effects numerically. 
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Figure 2. A known reference layer over the Cu-Ni sample with smooth potential at boundaries. The red 
lines demonstrate the ideal sample with non-smooth boundaries. 

3. Numerical example 

As a realistic example to study the stability of the reference method by using a smooth potential, we 
consider a sample with 20 nm thick Copper on 30 nm Nickel having a constant SLD of 6.52 x 10""^ nm~^ 
and 9.4 x 10""^ nm~^ as the unknown layer; and three separate reference layers with 20 nm Au, 15 nm 
Cr and 10 nm Co having a constant SLD value of 4.46,3.03 and 2.23 (x 10""^ nm"^), respectively. Figureji] 
shows this configuration for the gold layer as reference. Figure[3]also shows the reflectivity of the sample 
of figure[2]for three difference smoothness factors. The smearing of potential at boundaries is taken into 
account by using cr = 0, 0.2 and 0.5 nm. The effects of smoothness on output results are clear in the 
figure, particularly at large wave numbers, while the data for small wave numbers truly correspond with 
none-smooth data. 

Reflectivity of the unknown part of the sample with smooth varying SLD at the interface with a = 5 A 




0.4 0.6 

q (1/nm ) 



Figure 3. (Color online) Reflectivity from the sample of figure[2]for three cases: a = 0, 0.2 and 0.5 nm. 



13604-5 



S.S. Jahromi, S.F. Masoudi 



is illustrated in figure [4| (upper curve). Similar to figure |3] the effects of smoothness are more obvious 
at large wave numbers. However, consideration of smoothness at interface would cause some abrupt 
noises at certain values of wave numbers such as: q - 0.3 and 0.6 nm~^ The noises also appear for the 
real and imaginary data of reflection coefficient (figure [s] (a), (b)). These noises are due to the abrupt 
changes in some of the elements of the M matrix, equation l [T2) , in the reference method. The intensity 
and distribution of these noises along the different ranges of wave numbers are extremely relevant to 
the smoothness factor, thickness of the reference layers and the range of neutron wave numbers. In 
order to better understand the origin of these noises, we have plotted the M(3, 3) element of the M matrix 
in figure |4| (lower curve). The result shows that M(3,3) has some noises in q values, exactly as for the 
reflectivity data. Such severe noises in output results would cause a complete loss of reflection coefficient 
data in those q values. As we mentioned in section[2] the data of reflection coefficient in the whole range 
of wave numbers, even bellow the critical wave number, are needed to retrieve the SLD of the sample. 
Unfortunately, the noises are distributed over different values of wave numbers ranging from small to 
large q's,. This deficiency makes the data useless in retrieving the SLD of the sample. Here we are going 
to propose a method for removing or reducing these noises. 




- M(3.3} 



Figure 4. Reflectivity of the unknown sample of figure [2] which is reconstructed by using three reference 
layers: 20 nm gold, 15 nm Cr and 10 nm Co (Circled curve is for cr = 0.5 nm and solid line is for the ideal 
sample, a = 0). M(3,3) element of the M matrix, equation (l2j, is plotted for better understanding of the 
origin of reflectivity noises. 



The stability of the method to a great extent depends on the M matrix which is a function of the known 
part of the sample. The best result for the reflection coefficient of the unknown part of the sample comes 
out of the equation ( [12) in the case where the behavior of the M matrix is smooth. Our investigation 
shows that not every arbitrary reference layer can lead to a smooth result for the M matrix. Since we 
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can control the choice of the reference layers, we can run a simulation to find the best layers for which 
the M matrix behaves smoothly and without noise. By choosing the best reference layer, we can retrieve 
reliable results for the unknown part of the sample. Those reference layers which are well behaved can 
be used for experimental implementation of the reflectometry experiments. As we mentioned before, 
the noise strongly depends on the thickness of the reference layers, and the choice of a proper thickness 
for the reference films would enhance the accuracy of the output results. In our example, the numerical 
calculations show that the best result with the smallest noise level would be obtained for 15, 10 and 20 nm 
thicknesses for Au, Co and Cr, respectively. The purified curves of real and imaginary parts of reflection 
coefficient in this case are shown in figures |6] (a), (b). As it is clear in the figure, we have much better 
results with less noise in comparison with figure [s] 



T ' 1 ' 1 ' 1 ' r 




1 1 1 1 1 1 1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 

q(1/nm ) (^j 




non-smooth 
a=5 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 

q(1/nm ) (j^) 

Figure 5. The reflection coefficient of the unknown part of the sample of figure[2]for two different a = 0, 
5 A by three reference layers as; 20 nm Au, 15 nm Cr and 10 nm Co. (a) Real and (b), imaginary parts of 
the reflection coefficient. 

Thereupon, the choice of an appropriate thickness for the reference layers would completely remove 
the noise or displace it to large wave numbers, where it would cause no ambiguity in retrieving the SLD 
profile of the sample. The rest of the noises can be removed by extrapolation. By using these corrected 
data as input for some useful codes such as the one which is developed by Sacks 1 10 - 12 1 based on the 
GeTfand-Levitan integral equation 1 1 - 3 1, the scattering length density of the sample is retrieved. The data 
of reflection coefficient up to ^7 = 0.6 are sufficient to retrieve the SLD by Sacks code and we can neglect 
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the data of large wave numbers. Hence, after purification of data, any remaining noise in the range of 
large wave numbers can be neglected. 

The SLD profile corresponding to r^^iq) was retrieved by using the real and imaginary parts of the 
reflection coefficient at the presence of the smoothness. As it is shown in figurelv] the method of reference 
layers is stable at the presence of smoothness and the SLD corresponding to r^^q) is truly reconstructed. 
The circled curve depicts the retrieved SLD profile of the sample. The retrieved thicknesses corresponds 
well with the thickness of the sample. 



T 




I ' 1 ' 1 1 \ ' 1 ' 1 ' 1 
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Figure 6. The reflection coefficient of the unknown part of the sample of figure[2]for a = 0,5 A. The three 
reference layers are: 15 nm Au, 20 nm Cr and 10 nm Co. Due to the proper selection of reference layers, 
the noise is removed in comparison with the result of figure |4| (a) Real and (b), imaginary parts of the 
reflection coefficient. 

Consideration of smoothness of interfaces makes the output results more corresponding to experi- 
mental data. In order to verify the efficiency of the method, the reflectivity curves were plotted in figurejs] 
for a sample from |9 |, Al203(6.7 nm)/FeCo(20 nm)/GaAS, for either of smooth (green solid line) and non- 
smooth boundaries (red dashed line), in contrast to the experimental reflectivity curve. As it is shown 
in the figure, the reflectivity curve in the presence of smoothness of interfaces better corresponds with 
the experimental reflectivity curve while the simulated reflectivity for an ideal sample does not abso- 
lutely correspond with the experimental curve, particularly at large wave numbers. The results verify 
that consideration of the smoothness of interfaces is of highest importance. 
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Figure 7. The reconstructed SLD of the unknown sample shown in figure [2] surrounded by identical 
medium on both sides iQhh^, obtained by inversion of the extracted values of the reflection coefficient 
in the momentum range ^ ^7 ^ 0.6 nm~^ . The non-smooth profile is shown by thin solid line. The circled 
curve represents the SLD of the sample in the presence of smoothness. 



4. Conclusion 



In most of the simulation methods of neutron reflectometry, the SLD of adjacent layers is supposed 
to be discontinuous and sharp. Since in real samples there is observed some smearing at boundaries, 
consideration of smoothness would effect the output results and the retrieved data. 

In this paper by considering a certain type of Error function which denotes the smooth variation of 
SLD across the interface and implementing it into the reference method, the effects of smoothness on 
determining the phase and retrieving the SLD profile of the sample were investigated. Due to consider- 
ation of smoothness of interfaces, several abrupt noises existed in the output, such as reflectivity, real 
and imaginary parts of reflection amplitude. It was shown that these noises are due to the abrupt change 
in the elements of M matrix of equation ( [T2] l. We proposed to choose a proper thicknesses for the refer- 
ence layers, which would completely remove the noise or shifted it to the large wave numbers iq > 0.6), 
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Figure 8. (Color online) Experimental reflectivity curve in contrast to the simulated reflectivity at the pres- 
ence of smoothness of boundaries and the ideal sample. The green solid line demonstrates the reflectivity 
at the presence of smoothness and better corresponds with the experimental curve. 
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where the data of reflectivity are not needed in order to retrieve the SLD and we can easily neglect them. 
Consequently, the SLD of the sample was retrieved by using the purifled data of the reflection coefficient. 

In order to make sure that the method is highly reliable, we also compared the simulated reflectivity 
with the experimental curve for a sample. Astoundingly, the results show that consideration of smooth 
variation of SLD at the interface would make the measured reflectivity to better correspond with the 
experimental data. The results also certify that the reference method is stable in the presence of smooth 
interfaces. 
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Aoc/iiA>KeHH5i Bn/ii/iBy r/iaAKOCTi Mi>K(|>d30Bi/ix rpaHi/iL^b 

Hd CTIMKICTb 30HAyK)m/IX HdHOMdCLUTdGHI/IX TOHKI/IX n/IIBOK 

3d AonoMoroK) HeMTpoHHoY pe(|>/ieKTOMeTpM 

C.C. flxapoMi, CO. Macyfli 

OiBi/iHHi/iM (t)aKy;ibTeT, TexHO/ioriMHi/iM yHiBepci/iTei iM. K.H. Tyci, a/c 15875-4416, TerepaH, IpaH 

Bi/ibLuicTb MBTOAiB pec|);ieKTOMeTpiT, 1140 BHKopncTOByK)Tbc?i fl,m BH3HaMeHH?i c|)a3n KOMn/ieKCHoro Koe(t)ii4ieH- 
Ty Bi^GnBaHMfi laKi 9\k era/iOHHUM mbtoa i 3MiHa npn/ier/ioro cepeAOBHii^a rpyHTywTbCJi Ha posB'jiBKy piBHfi- 
HH9\ LUpeAiHrepa 3 BUKopucraHH^iM pospuBHoro i cxoAHHKonoAi6Horo onrnHHoro noreHi^ia/iy po3ciK)BaHHfi. 
ripoTe, n'\/\ Mac npoi^ecy HanopoLuyBaHH^i fl,m ni/^roTOBKn pea/ibHoro spasKa ^Ba cyciAHi Luapn 3MiLuyK)Tbc?i i 
Mi>K(t)a30Ba rpaHi/ii4?i Moxe He 6yTn pospuBHOK) i HiiKOK). Po3MHBaHH?i cyciAHix LuapiB npn Mi>K(t)a30BiM rpa- 
Hi/ii4i (MaAKicTb iHTepc|)eMcy), Moxe Main Bnn\AB na B\fl,6\ABaHH9\, (^asy Koe(t)ii4ieHTa BiA6i/iBaHHJi i nepe6yA0By 
rycTMHM AOBXMHi/i po3ciK)BaHHfi 3pa3Ka. B i4iM cTaTTi MM /\ocj]\/\]Aj]]A cTiMKicTb eTa/ioHHoro MeTOAy y npi/icyTHocTi 
MaAKi/ix Mix(|)a30BHx rpaHi/ii4b. F/iaAKicTb Mi>Kct)a30Bi/ix rpahMi^b po3rn9\fl,aeThC9\, BMKopi/iCTOByKDHH HenepepB- 
Hy (t)yHKi4iK) noreHi^ia/iy po3ciK)BaHH?i. Mn laKox 3anponoHyBa/ii/i MeroAA-^^ 0Tpi/iMaHH?i HaM6i/ibLu HaAiwHoro 
pe3y;ibTaTy, fiKnw Bi^HOB/iKDe ryciMHy aobxhhi/i po3ciK)BaHHfi 3pa3Ka. 

K/iiOMOBi c/iOBa: nempoHHa pecp/ieKTOMerpm, tohki' nniBKi/i, r/ia/^Ki/iM noTeHU,ian, era/iOHHi/iyi Mewfl, 
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